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ABSTRACT 

The cosmological expansion cannot produce the reported correlations of 
the gamma-ray burst timescale and spectral energy with peak flux if the burst 
model reproduces the BATSE 3B peak-flux distribution for a non-evolving burst 
source density. The required ratios of time dilation and redshift factors are only 
produced by monoluminous models at peak fluxes below the BATSE threshold, 
and they are never produced by power-law luminosity models. Monoluminous 
models only produce acceptable fits to the peak-flux distribution for very 
specific combinations of the spectral and cosmological parameters. The redshift 
of gamma-ray bursts at the BATSE threshold is z ~ 1.5. Power law luminosity 
distribution models oc produce acceptables fits to the data for most values 
of the spectral parameters when (3 < 1.6. In this model, gamma-ray bursts of 
a given peak flux have a distribution of redshifts, with a maximum redshift of 
^ 3 for peak fluxes near the BATSE threshold, and with an average redshift 
of < 1 for all values of peak flux. This qualitative behavior occurs whenever 
the luminosity distribution determines the shape of the peak-flux distribution, 
regardless of whether source density evolution occurs. The reported correlations 
of the burst timescale and the spectral energy with peak flux are systematically 
1 standard deviation above the monoluminous model and 1.5 to 2 standard 
deviation above the power-law luminosity model. These results suggest that 
an intrinsic correlation of burst timescale and spectral energy with luminosity 
is present. Studies of the peak-flux distribution for bursts selected by E peak or 
hardness ratio provide a test for this intrinsic correlation. 
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1. The Indistinctness of Cosmological Signatures 

Gamma-ray bursts are widely believed to be cosmological, but no direct measurement 
of distance exists to confirm this belief. No gamma-ray burst has lines that can be associated 
with atomic or nuclear processes, so a direct calculation of their redshifts is not possible. 
With only one recent exception (Groot et al. 1997), no gamma-ray burst has an observed 
counterpart, so an indirect calculation of their redshifts is also not possible. Because of 
these difficulties, a number of researchers have attempted to derive average redshifts by 
studying the properties of gamma-ray burst catalogs. The simplest of these studies fit the 
peak-flux distribution to cosmological models under various assumptions (e.g. Fenimore 
et al. 1993; Wickramasinghe et al. 1993; Emslie & Horack 1994; Cohen & Piran 1995; 
Pendleton et al. 1996; Petrosian & Lee 1997). More involved studies attempt to measure a 
correlation of gamma-ray burst timescale with photon flux under the assumption that such 
a correlation is evidence of a cosmological time dilation (Norris et al. 1994; Davis et al. 
1994; Mitrofanov et al. 1994, 1996 Lee & Petrosian 1996). Related studies attempt to 
measure a correlation of the characteristic photon energy with flux, which may indicate a 
cosmological redshift (Nemiroff et al. 1994; Mallozzi et al. 1995). 

The difficulty with these approaches is that the interpretation of the results is not 
unique. The observed flux distribution simply shows that gamma-ray bursts are not 
homogeneously distributed in space, which is unique neither to a cosmological origin nor 
to a particular model of cosmological expansion. Galactic corona models of gamma-ray 
bursts easily reproduce the observations (Brainerd 1992; Ulmer & Wijers 1995; Ulmer 
et al. 1995; Hakkila et al. 1995), and any cosmological model can be made to conform to 
the observations with the proper choice of density evolution (Rutledge et al. 1995; Horack 
et al. 1995; Reichart & Meszaros 1997), or luminosity distribution (Meszaros & Meszaros 
1995; Hakkila et al. 1996), as is evident from the numerous combinations of assumptions 
made by the numerous papers examining this topic. The effects sought in correlation 
studies can be explained by an intrinsic correlation of the parameter of interest with the 
luminosity (Brainerd 1994a). The beaming of relativistic jets (Brainerd 1994a; Yi 1994) 
and the conservation of total energy in monoenergetic sources (Wijers & Paczyhski 1994) 
are examples of mechanisms that produce such a correlation. 

The question addressed in this article is whether a strong correlation of the burst time 
scale or the characteristic spectral energy with peak flux can have a purely cosmological 
origin when the gamma-ray burst peak-flux distribution is determined by the luminosity 
distribution rather than by the spatial distribution. In this study, I use a luminosity 
function that is a power law with an upper limit, and, for comparison, a monoluminous 
distribution. The spectral model I use (§2) is the Compton attenuation model (Brainerd 
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1994b), which accurately represents the observed spectra. The cosmological expansion 
with a (1 + zY comoving source density (/i < 9/2) defines a limiting power-law index 
for the peak flux distribution that is > —5/2 as the peak flux goes to zero; a luminosity 
distribution can drive this index closer to —5/2 (§3 and §4). When the luminosity function 
determines the asymptotic behavior of the peak flux distribution, the density evolution 
of the gamma-ray burst sources has only a weak effect on the average time-dilation and 
redshift (§3). I find (§5) that a broad luminosity distribution in a Friedmann cosmology 
with a constant comoving source density produces peak-flux distributions that fit well the 
distribution observed by the Burst and Transient Source Experiment (BATSE) on the 
Compton Gamma-Ray Observatory. At a given flux, the gamma-ray bursts have a redshift 
distribution, with the average redshift considerably smaller than the maximum redshift (§4). 
The maximum redshift a burst can have at the BATSE threshold is ^ 3. In this model, 
the cosmological expansion cannot explain the reported correlations (§6), so an intrinsic 
correlation must exist. The flux distribution for gamma-ray bursts selected by hardness 
ratio presented in recent studies is a measure of this intrinsic correlation (§7). These studies 
offer the opportunity to separate the cosmological effects from the intrinsic correlations. 

2. Spectral Parameterization 

Gamma-ray burst spectra are not power-law spectra; they are hard spectra that 
gradually roll over in the 50 keV to 300 keV energy band to a power law, with the power-law 
behavior appearing above 1 MeV (Metzger et al. 1974; Hueter 1984; Matz et al.1985). The 
choice of spectral model influences the shape of the burst peak-flux distribution (Brainerd 
1994b; Mallozzi et al. 1996). A physical spectral model that accurately describes gamma-ray 
burst spectra is a Compton-attenuated power-law spectrum (Brainerd 1994b; Brainerd 
et al. 1996a,b). Because this spectrum is a testable gamma-ray burst theory, it is of interest 
beyond its descriptive value. 

This spectrum is defined by four parameters in the gamma-ray energy band: a 
Thomson optical depth tt, a parameter \1> = (5 + 2) jr T < 0, where 5 is the power-law 
index of the spectrum as v — > oo, a redshift z, and an overall normalization. In the 
discussion that follows, the spectral parameter z is set to 0, so that the cosmological redshift 
enters explicitly through equation (7). The overall normalization is folded into the free 
parameter L, which is a luminosity that is defined in equation (6). The two remaining 
parameters describe the shape of the spectrum. The parameter \1/ sets the peak energy 
E p of vF v in the cosmological comoving frame, and the parameter tt defines the width 
of vF v about E p . The energy E p only exists for > —0.1972, taking on a minimum 
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Fig. 1 — Compton attenuated 
spectrum. For all three curves, 
the spectral parameter 5 = —4. 
From top to bottom, the curves 
have tt = 10, 15, and 20, 
corresponding to * = —0.20, 
—0.13, and —0.10 respectively. 



Energy (keV) 

value of 231 keV at * = -0.1972. For < < 0.1972, the value of E p is given by 
^ = — log (2E p /m e c 2 ) m e c 2 /E p , so a wide range of ^ values produces a narrow range 
of -Ep values. The width in logE of the vF v curve about E p is proportional to 1/y/vr- 
General values of these parameters that produce good fits to the data are 10 <^ tt <^ 30 
and —0.16 < ^ < —0.12 (Brainerd et al. 1996b). Three examples are given in Figure 1. In 
this figure, 5 = -4 and r T = 10, 15, and 20, which gives ^ = -0.20, -0.133, and -0.10 
respectively; the uppermost curve does not have an E p since ^ < —0.1972. The inflection 
point at 231 keV is apparent in these spectra. The x-ray upturn predicted by this spectrum 
has been observed (Preece et al. 1996; Brainerd et al. 1996a,b). In this study, only the 
spectrum above 50 keV is used. 



3. The General Gamma-Ray Burst Flux Distribution 

The equation for the flux distribution of cosmological gamma-ray bursts is well known. 
The general equations for the number of gamma-ray bursts per unit proper time r per unit 
proper volume V per unit luminosity L per unit spectral parameter space p is 

dN 



dr dV dL d n p 



p(T,r,L,p) = R(t) n (r)$(L,p] 
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where the second equality expresses the assumption that the characteristics of a gamma-ray 
burst are independent of the age of the universe. The function n (r) is the number of 
gamma-ray bursts per unit proper time per unit comoving volume. The function $ (L,p) 
is the fraction of gamma-ray bursts per unit luminosity per unit spectral parameter space; 
it is normalized to unity. The proper time of the gamma-ray burst source is related to the 
observer's time by 

dt = (1 + z) dr, (2) 
and the volume element is related to the coordinate elements through 

2 

dV = R{rf ^drdtl, (3) 

v /l-(2g -l)r 2 

where the nonstandard coordinate radius r is defined in terms of the Robert son- Walker 
coordinate radius r\ as r = riH R /c. The parameter R(r) is c/H Ro times the expansion 
parameter for the Robertson- Walker metric, and the factor Jl — (2q — l)r 2 in equation 



(3) is the familiar ^/l — kr\ from the Robertson- Walker metric. The coordinate radius is 
related to the redshift by 



q z + (go - 1) (V 1 + 2( loz - l) 2z ( 1 + z + VI + 2%; 



* 2 (1 + ^) (l + *)(l + Vl + 2gb*)' 



(4) 



From the formal derivation of equation (4) (Weinberg 1972, eq. [15.3.23]), the first derivative 
of r with z is 

dr_ = ^l-(2g -l)r 2 

dz V 1 + 2?o-2 (1 + z) ' U 

The gamma-ray number flux F in an energy band spanning the frequencies V\ to v 2 is 
related to the luminosity by 



LH, 



F = ^F(z,p), (6) 
<r 

where the length scale Ro is folded into the luminosity L, and the function T which 
is defined as 

F(z,p)=r- 2 I[u(l + z),p] du, (7) 

J V\ 

is the photon number flux normalized by the overall luminosity. The function I (u) is the 
photon number flux per unit frequency. The purpose of this particular parameterization is 
to separate the luminosity L from the spectral shape parameters p. 
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With equations (2) through (7), one can transform equation (1) into a distribution 
over observable parameters: 

dN ( F \ r 2 

This equation differs from the equation for galaxy distributions by a factor of 1/ (1 + z), 
which arises from the time dilation of the burst rate. 

The distribution of bursts with flux is given by the integration of equation (8) over all z. 
For this integral to be bounded, $ (L,p) must have limits on its asymptotic behavior as 
L — > and L — > oo. If we assume that $ (L,p) oc L~& 1 for L <C L , where L characterizes 
the peak of the function L 5//2 $ (L,p), then the integral of equation (8) over z converges as 
z — > when $ < 5/2. If $ oc L _/3/t for L 3> L , the integral converges as z — > oo 

when /3ft, > 5/2. While the second limit is derived under the assumption that 2; < 1 when 
L ^> Lo, relaxing this assumption lowers the value on the right side of the inequality, so the 
condition holds for all parameter space. 

3.1. Asymptotic Behavior 

The flux appears in equation (8) only in the argument of the luminosity distribution 
function. The asymptotic behavior with flux of the integral of equation (8) over z is 
therefore determined by the behavior of $ (L,p) at the values of z producing the greatest 
contribution to the integral. If this occurs for values of z that give L ~ Lq, then the 
integration variable z can be replaced by L, and the expansion of factors involving z 
determine the asymptotic behavior with F. On the other hand, if the integrand peaks at 
values of z such that L <C L or L ^> L Q , then it is the asymptotic behavior of $ (L,p) that 
determines the asymptotic behavior of the flux distribution with F. 

When z 1 for L ^> Lq, the flux is oc L/z 2 , and equation (8) is oc $ (L,p)z 4 , 
assuming no density evolution at z <S 1. The integrand peaks at L , so its contribution 
to the integral above z ~ 1 is negligible. Integrating over z and changing the variable of 
integration to L then gives the flux distribution, which is proportional to F~ 5 ^ 2 . This is the 
origin of the —3/2 slope of the logiV-logF curve of a homogeneous distribution. 

Of more interest is the behavior of equation (8) when z ^> 1 for L ~ L . To discuss the 
asymptotic behavior of the integral in this regime, I let the photon spectrum be a power 
law, so that l v oc u 6 , with 5 < —2. This is justified by the observed power-law behavior 
above 1 MeV. I also assume that $ (L,p) oc L^ 13 when L > L or L C L , with (3 = (3 h for 
the first condition and P — Pi for the second. Finally, I set n (r) oc (1 + zY to investigate 
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the effects of density evolution. A necessary and sufficient condition for the integral over 
equation (8) to be bounded under these assumptions is that (5\ < 5/2, which guarantees 
the convergence of the integral as z — > 0, and fa > 1 + (3/2 — /x) /S, which guarantees 
convergence as z — > oo. Because fa > 5/2 is already required for convergence in the regime 
of zo <C 1, this last criterion is satisfied when /x < 3 (1 — 5)/2. Because 5 < —2, the right 
side of this inequality is > 9/2, so this second condition on fa is only important for very 
strong evolution. I assume that this limit on /x holds in the discussion below. 

When z ^> 1 for L ps Lo, the integral can be divided into four parts that are bounded 
by z = 1, z = 1/qo > 1, and L (z) = Lq. When the integral peaks at a value zq that is either 
of these first two values, the factor of F in the argument of the luminosity function can be 
factored out of the integral. The distribution in flux is then proportional to F^ 13 . There 
remains two questions to answer: when is z at L (zo) ~ L , and what is the dependence 
on F of the flux distribution in this case? 

Two regimes exist that affect the dependence of the burst distribution when L (z ) ~ L : 
1 < Zq < l/q and z > 1/qo. In the first regime, equation (8) has the asymptotic integral 
form of 

dN f°° ^ ( F J\ _ , v_i „ , 

,p) F{z) 1 dz. (9) 



dtdQdFd n p Jo \F (z 
Because T {£) oc z~ 2 ~ 5 in this regime (eq. [7] and [4]), this equation becomes 

dN „ i i+M 

oc i?- 1 -— . (10) 



dt dtt dF d n p 



Necessary and sufficient conditions for L(z ) L are f3i < 1 + |±| and fa > 1 + . 
Because I am assuming /x < 9/2 and 5 < —2, the strongest condition on fa is /3ft > 19/11, 
which is weaker than the earlier condition of fa > 5/2. As a consequence, one never can 
have z ~ 1/qo- The only possibility is to violate the condition on /3j, in which case z ~ 1 
in the 1 < z < l/q regime. Because of this, the general asymptotic equations for the burst 
flux distribution is 

dN 

ocF~ ai (11) 

dtdndFd n p ' { } 

where 

ai = max (ft, 1 + • (12) 

In the second regime, where z > 1/qo, the asymptotic form of the integral over 
equation (8) is 

dN r°° ( F \ ib 

n c * hF7T>P Fiz)- 1 zT* dz. (13) 



dt dQ dF d n p 
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In this regime, J- (z) oc z , so the dependence of the flux distribution on F is 

oc F 1+ s . (14) 



di dfi dF d n p 




Sufficient (but not necessary) conditions for L (z ) L are $ < 1 + Pi < 1 — M ~|/ 2 , 
and // > (6 — 55) /4. The last condition guarantees that the integration does not have 
secondary maximum at z — 1. Because I assume 5 < —2 and /i < 9/2, this last condition 
is always satisfied. As a consequence, the second condition on (3i is the stronger, and the 
general asymptotic solution can be written as 

dN 

oc F~ a2 (15) 

dt dn dF d n P ' { } 

where 

/ „ - ii \ 

(16) 

Several points are immediately clear from these equations. First, the behavior of 
the luminosity distribution function below Lq influences the asymptotic behavior of the 
flux distribution, while the behavior above Lq is unimportant. Recall that Lq is defined 
as the maximum of L 5//2 $ (L,p). Second, in the absence of strong source evolution, a 
large asymptotic power-law index at small F requires a luminosity distribution with a 
power-law index of equal value, as noted by Meszaros & Meszaros (1995). Because 5 < —2, 
the minimum value of cci, which is given by the second term in equation (12), is < 5/4 
for \x = 0; from the second term of equation (16), the minimum value of «2 is < 1/4 for 
\i = 0. Third, because the luminosity function can formally diverge and still give a bounded 
integral as long as it falls below a L~ 5 / 2 power law as L — > 0, the presence of a change in 
the luminosity distribution function at sufficiently low luminosity — which is required by 
the boundedness of / $ (L) dL — has no effect on the observed flux distribution; formally, 
this must be at L F min c 2 / HqJ 7 (1), where F min is the minimum observable burst flux 
(Meszaros & Meszaros 1995). Finally, evolution is unimportant in the asymptotic behavior 
unless the number of bursts increases dramatically as 1 + z increases; the second index 
in equation (12) increases by < /i/4, and the second index in equation (16) increases by 
< /i/2, so evolution cannot dominate the luminosity distribution function unless ji ^ 3. 

Equations (12) and (16) imply different behaviors for q = 0.5 and q <^ 0.5. For 
q = 0.5, the only relevant limit is that given by equation (16). When z 1, the integrand 
rises with z until the peak in $ (L,p) is reached. When zq 1, the integrand rises until 
z w 1; then, if ft > 1 — (jjl — §) /5, the integrand goes to zero more rapidly than z^ 1 
as — > oo, but if the inequality is violated, the integrand either rises or goes to zero less 
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rapidly than z^ 1 . In contrast, when q <C 0.5, there are three regimes to consider. As in 
the previous case, when z < 1, the integrand rises to the maximum of $(L,p), but for 
1 < z < 1/qo, equation (12) applies, and the integrand goes to zero more rapidly than z~ l 
when Pi > 1 + (1 + /i) I (2 — 5), and rises or goes to zero less rapidly than z~ l when the 
inequality does not hold. When zq > 1/qo, equation (16) holds, and the behavior noted for 
go = 0.5 holds for q <C 0.5. The consequences of this behavior are that the minimum value 
of a in an F~ a asymptotic description of the flux distribution is strongly dependent on q 
and that, for q <C 0.5, the distribution in z for bursts with a given F is broad and flat over 
1< z < l/q when 1 + (1 + //) / (2 - 5) > ft > 1 - (fi - §) /5. 



4. Incorporating a Luminosity Distribution 

The gamma-ray burst flux distribution is derived for two luminosity distributions: a 
monoluminous distribution function, and a power-law luminosity distribution with a high 
luminosity cutoff. The first is widely used, and it is given here for comparison to the power 
law distribution. 



4.1. Monoluminosity Distribution 

The equation for a monoluminous gamma-ray burst distribution does not have an 
integral over z. It is therefore relatively easy to evaluate and fit to the observed flux 
distribution. From equations (1) through (6), 



dN 4r 2 
- iV 0Tr- — ^2" 



dF (l + zf r^l + 2q z [ dT{z)jdz 



dT [z] jdz\ 



(17) 



In this equation, tq is the value of r from equation (4) for z — 1. The normalization Nq is 
the number of bursts per unit F at z — 1. F is related to z through 



'JF(l)' vxu; 

where T (z) is given by equation (7). 

The distributions produced by equation (17) are plotted as dotted lines in Figure 2 for 
F = 1, iVo = 1, ^ = —0.15, and t t = 15 for both q = 0.5 and 0.1. The curves for z given 
by equation (18) as a function of F are plotted in Figure 4 for these parameters. The curves 
in Figure 2 go to the correct limiting power-law indices: for q = 0.5, it is given by equation 
(16) as 11/17 (w 0.65), while for q = 0.1, it is given by equation (12) as 29/25 (w 1.16). 
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Flux Flux 

Fig. 2 — Peak-flux distributions. All curves are for spectral parameters 
^ = —0.15 and r T = 15, which are defined in §2. The figures are for q = 0.5 
(left) and 0.1 (right). The solid curves for a power-law luminosity distributions 
with (top left to bottom left) (3 = 1.8, 1.4, 1.0, and 0.6. The solid curve 
is the monoluminosity distribution for the given spectral and cosmological 
parameters. The curves are normalized by iVo = 1 and Fq = 1 so that 
dN/dF = 1 and z max = 1 at F — 1. The flux is measured in the 50keV 
to 300 keV energy band. 



4.2. Cutoff Power-Law Luminosity Distribution 



A simple distribution function used frequently in gamma-ray burst studies is a 
power-law distribution with one or two luminosity cutoffs. In this study, I use 



$ (L) = < 



(£) _/5 , forL<L ; 
, , otherwise. 
Placing this into equation (8) and integrating over z gives 



(19) 



^l = n U-) / T^ x {z) J dz, (20) 

dF °\F J Jo V ; (1 + zf y/l + 2q z ' 1 ' 

where z max is defined for a given flux F by equation (18) with z replaced by z max . The 
normalizing flux F is the flux at which z max = 1- The normalization parameter n is 
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defined as 

1 2 

n = N / [ T^ 1 {z)——^-==dz, (21) 
Jo (1 + 2) V 1 + 2 <?0<2 

where N is the number of gamma-ray bursts per unit flux at F = F . 

Flux distributions from equation (20) with F = 1 and N = 1 are plotted in Figure 2 
as solid lines for (3 = 0.6, 1.0, 1.4, and 1.8 for both go — 0.5 and 0.1. The asymptotic 
behavior discussed in §3.1 is clearly demonstrated in these figures. For go — 0.5, equation 
(16) states that the asymptotic power-law index is given by —/3 when (3 > 11/17 ~ 0.65. 
Three of the four solid curves in this plot satisfy this criterion, and all three are close to 
their asymptotic values at low flux. The fourth, which has (3 = 0.6, is slowly going to 
the asymptotic value of —(3 = —11/17. For g = 0.1, the limit on (3 is given by equation 
(12), because z max > 6 for the fluxes in the plot, and equation (16) is not applicable until 
z 3> 10. In this case, the asymptotic limit is given by —f3 when (3 > 29/25 = 1.16. In 
Figure 2b, the upper two solid curves go to a power-law index of —(3 as F — > 0, and the 
lower two solid curves approach the monoluminous curve. 



4.3. Average Redshift 

An important difference between a monoluminous distribution and a broad luminosity 
distribution is that in the former the redshift and the flux have a one-to-one correspondence, 
while in the latter there is a distribution of redshift values for each flux value. From 
equation (20), one sees that the integrand is independent of F, and that the dependence 
on flux of the integral arises from the dependence of z max on flux in equation (18). The 
shape of the redshift distribution is therefore independent of flux, and only the upper 
limit on z changes with flux. When the integrand rises monotonically with z for all z, the 
median value of z is a strong function of z max , but when the integrand falls sufficiently 
rapidly above some value of z, the median value of z goes to a constant as z rnax goes to 
infinity. From the discussion of the asymptotic behavior of equation (8), we know that the 
first instance occurs when L(z ) L , in which case the dependence of the burst number 
distribution on F is determined by the second index in either equation (12) or equation 
(16). In this case, the asymptotic behavior of the flux distribution is independent of the 
luminosity distribution function. The second instance arises when zq ~ 1 or zq ~ 1/go, i n 
which case the asymptotic behavior of the flux distribution is determined by the luminosity 
distribution. 

Figure 3 shows the z distribution function for (3 = 0.6, 1.0, 1.4, and 1.8. In the cases 
where (3 is large, the distribution peaks below z = 1 and falls rapidly as z increases, so 
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Fig. 3 — Distribution of z for a power-law luminosity distribution. On each 
figure, the distribution curves for power-law indexes (from upper left to lower 
left) of (3 — 1.8, 1.4, 1.0, and 0.6 are plotted as functions of cosmological 
redshift z. The curves are normalized to unity at z — 1. The figures are for 
q = 0.5 (left) and 0.1 (right). The spectral parameters are \& = —0.15 and 
tt = 15. 



the average value of z is a weak function of z max when z max ^> 1. As (3 decreases, the 
peak in the distribution function moves to higher z, until it disappears altogether, and 
the distribution function rises monotonically with z, so that the average value of z goes to 
Zmax- Comparing Figure 3a to Figure 3b, one sees the behavior described in §3.1. When 
qo = 0.5, the integrand, and therefore the distribution in z, falls faster than z^ 1 as z — > oo 
when (3 > 11/17. The three curves that satisfy this criterion exhibit the expected behavior. 
The curve that violates the inequality falls slowly as z rises. When q = 0.1, the criterion is 
(3 > 29/25, which only two of the curves satisfy. Of the two curves violating the criterion, 
the (3 = 1 curve falls slowly, while the (3 = 0.6 curve continues to rise. When z > 10, which 
is not shown in the figure, both of these curves fall, since the (3 > 11/17 inequality from 
equation (12) become the governing criterion. 

Tests for the effects of the cosmological expansion rely on the average properties of 
gamma-ray bursts. Because the bursts of a given peak flux have a distribution of redshifts, 
a distinction must be made between the value of (z) derived from time-dilation and the 
value derived from redshift. For the first, one is averaging 1 + z, while for the second, one is 
averaging 1/ (1 + z). The consequence is that the average value of z derived in the second 
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Flux Flux 

Fig. 4 — Average values of z as a function of flux. The solid line gives the 
relationship between z and F from equation (18) with F — 1. The dotted 
lines give (z), while the dashed lines give l/^l + z) -1 ^ — 1. The different 
lines are for, from bottom to top, (3 = 1.8, 1.4, 1.0, and 0.6. The figures are 
for g = 0.5 (left) and 0.1 (right). The spectral parameters are \I/ = —0.15 and 
tt = 15. 

case is smaller than in the first. This is shown in Figure 4, where (z) is plotted as dotted 
lines, and 1/ ^(1 + z) -1 ^ — 1 is plotted as dashed lines for the values of (3 given in Figure 2. 
The curves go to a constant as z — > oo when they satisfy the criterion on (3 that makes —f3 
the asymptotic power-law index of the flux distribution. 



5. Fitting the BATSE Flux Distribution 

I use x 2 minimization to fit the model flux distributions given in §4 to the flux 
distribution for a subset of gamma-ray bursts from the BATSE 3B catalog. This catalog 
gives the peak fluxes in the energy range of 50 keV to 300 keV for 869 gamma-ray bursts 
on three different time scales: 64 ms, 256 ms, and 1024 ms. The subset I use comprises the 
bursts that are not overwrites of a preceding burst, that have a T90 duration > 2s (where 
T90 is the time over which 90 % of the burst energy is emitted), and that have a peak flux 
on the 256 ms timescale of 0.7cm -2 s -1 < F < 110.9 cm~ 2 s _1 . There are 396 gamma-ray 
bursts that meet these criteria. 
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The upper threshold is for computational convenience, and it has no effect on the 
model fits. The lower threshold is to avoid the effects of the variable detector threshold 
and the detector geometry. The BATSE instrument triggers on a gamma-ray burst if the 
count rate exceeds a threshold count rate in at least two detectors. This count rate is set 
by integrating the background for 17 s intervals and calculating from this the standard 
deviation for counts on each of the three timescales given above. The trigger threshold is 
a multiple of this standard deviation; usually 5.5 standard deviations is chosen. Because 
the background fluctuates by a factor of 2 over the 90 minute orbit, the trigger threshold is 
variable. Because the trigger is in counts rather than in flux, because the detector geometry 
is an octahedron, and because Earth scatters burst radiation to the experiment, the flux 
threshold is angle dependent. These various effects cause an undercounting of gamma-ray 
bursts of low flux. The choice of 0.7cm _2 s _1 on the 256 ms timescale as a lower limit 
overcomes these difficulties (Pendleton et al. 1996). 

For 90 minutes after a trigger, the threshold is raised so that if a particularly bright 
burst occurs, it is observed by the instrument. Such bursts, which are called overwrites, 
must be dropped from the sample, because they represent an overcounting of bright bursts 
relative to dim bursts. 

The choice of timescales affects the choice of cosmological model. The reason is that 
if the burst pulse carrying the peak flux is fully resolved in time, then the cosmological 
redshift comes into the peak flux through the Lorentz boost of the flux, but if the pulse 
is much shorter than the instrument timescale, all photons in the pulse arrive inside that 
timescale regardless of time dilation. The time dilation of the emission therefore drops 
out of the problem, and the peak-flux has one more factor of (1 + z) than it does for the 
time-resolved pulse. I desire a timescale that resolves the time-behavior of most gamma-ray 
bursts without introducing large errors in the flux measurement from counting statistics. I 
compare in Figure 5 the peak fluxes of the 396 bursts in my sample. If two time scales each 
resolve the time-dependent variation of the flux, the peak flux measured on each should be 
equal within counting statistics. In Figure 5a, each burst is plotted on the 256ms-1024ms 
timescale plane, and in Figure 5b, each is plotted on the 64ms-256ms plane. The upper 
diagonal line represents equal fluxes on both timescales, which occurs when a pulse is fully 
time resolved. The lower diagonal line represents the short timescale having a peak flux 
that is 4 times the peak flux for the long timescale, which is what occurs when the pulse 
is much shorter than the short timescale. One sees that the clustering along the line of 
equality is stronger for Figure 5b than for Figure 5a, showing that the 256 ms timescale 
resolves most burst while the 1024 ms timescale leaves a large fraction of bursts unresolved. 
The percentage of bursts in Figure 5a that have a 256 ms peak flux that is more than a 
factor of 2 larger than the 1024 ms peak flux is 2.6%, while the percentage in Figure 5b that 
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256 ms Peak Flux (counts cm~ 2 s~ 1 ) 64 ms Peak Flux (counts cm~ 2 s~') 



Fig. 5 — Peak flux scatter plot. The gamma-ray bursts from the BATSE 3B 
catalog that have a duration measure T90 > 2 s are plotted by their (left) 256 
and 1024 ms peak fluxes and their (right) 64 and 256 ms peak fluxes in the 
50keV to 300 keV energy band. The upper diagonal line on each plot gives 
flux equality, while the lower diagonal on each gives a ratio of 1 to 4 in the 
two peak fluxes. 

have a 64 ms peak flux that is more than a factor of 2 larger than the 256 ms peak flux is 
0.7%. For a factor of y/2, these percentages are 19.8% and 10.2% respectively. This is why 
the 256 ms timescale is used. 

Some bursts have durations that are of order the sampling timescale. For these bursts, 
the time structure cannot be resolved, so a lower limit on the gamma-ray burst duration is 
warranted. The duration distribution is bimodal, with a local minimum at 2 s. This value 
is therefore a reasonable lower limit on T90 in the sample. 

The burst distribution is summed into 22 bins that uniformly spanning log F between 
0.7cm _2 s~ 1 and 110.9 cm -2 s" 1 . Then, going from the highest flux bin to the lowest flux 
bin, each bin with fewer than 20 bursts is added to the bin immediately lower in flux. The 
purpose of this is to make gaussian statistics the correct description of the number of bursts 
in each bin. This lowers the number of bins to 11. 
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Fig. 6 — Reduced x 2 f° r 
fits of monoluminous peak-flux 
distributions to the BATSE 
data. The x-axis is the spectral 
parameter tt, and the various 
sets of points are for ^ = 
-0.10 (+), -0.12 (*), -0.14 
(O), -0.16 (A), and -0.18 (□). 
In the upper figure, q = 0.5, 
and in the lower figure, go = 0.1. 



10 12 14 16 18 20 22 24 26 21 



-17- 



TABLE 1 

MONOLUMINOSITY, q = 0.5 







Reduced x 2b 


F C 


<?F d 


1 + 

1 + Z10 


— u. 1U 


1U 


1 .u 


1 Q 
1 .O 


u.o 


1 7fi 
1 . i t> 


-0.10 


12 


i.i 


1.3 


0.3 


1.80 


-0.10 


14 


1.3 


1.3 


0.3 


1.85 


-0.10 


16 


1.7 


1.4 


0.3 


1.90 


-0.10 


18 


2.2 


1.4 


0.3 


1.95 


-0.10 


20 


2.9 


1.4 


0.3 


2.01 


-0.10 


22 


3.8 


1.4 


0.3 


2.08 


-0.10 


24 


4.9 


1.4 


0.3 


2.15 


— U. 1 Z 


i n 
±u 


O.o 


^ A 


n q 
u.y 


z.uo 


-0.12 


12 


3.4 


1.8 


0.5 


1.88 


-0.12 


14 


1.4 


2.1 


0.4 


1.99 


-0.12 


16 


1.3 


2.2 


0.4 


2.05 


-0.12 


18 


1.2 


2.2 


0.4 


2.11 


-0.12 


20 


1.3 


2.2 


0.4 


2.18 


-0.12 


22 


1.4 


2.2 


0.4 


2.26 


-0.12 


24 


1.6 


2.2 


0.3 


2.35 


n i a 

— U. 


i n 
±u 


1 7 


1 7 


n a 


1 77 


-0.14 


12 


1.5 


1.7 


0.4 


1.79 


-0.14 


14 


1.4 


1.7 


0.4 


1.82 


-0.14 


16 


1.2 


1.7 


0.4 


1.85 


-0.14 


18 


1.1 


1.8 


0.4 


1.88 


-0.14 


20 


1.1 


1.8 


0.4 


1.91 


-0.14 


22 


1.1 


1.8 


0.4 


1.94 


-0.14 


24 


1.1 


1.8 


0.4 


1.97 


— U. ID 


1U 




9 Q 


U. t 


1 CI 


-0.16 


12 


3.2 


1.7 


0.5 


1.75 


-0.16 


14 


2.2 


1.9 


0.4 


1.78 


-0.16 


16 


2.1 


1.9 


0.4 


1.80 


-0.16 


18 


1.9 


1.9 


0.4 


1.81 


-0.16 


20 


1.8 


1.9 


0.4 


1.83 


-0.16 


22 


1.7 


1.9 


0.4 


1.84 


-0.16 


24 


1.6 


1.9 


0.4 


1.86 


-0.18 


10 


1.5 


1.1 


0.3 


1.61 


-0.18 


12 


1.5 


1.1 


0.3 


1.61 


-0.18 


14 


1.5 


1.1 


0.3 


1.61 


-0.18 


16 


1.4 


1.1 


0.3 


1.61 


-0.18 


18 


1.4 


1.1 


0.3 


1.61 


-0.18 


20 


1.4 


1.1 


0.3 


1.61 


-0.18 


22 


1.3 


1.1 


0.3 


1.62 


-0.18 


24 


1.3 


1.1 


0.3 


1.62 



a Spectral parameters that are defined in §3 and held 
constant in the x 2 minimization. 
b For 9 degrees of freedom. 

c The best fit value of the peak flux normaliza- 
tion in the 50keV to 300keV energy band in units of 
photonscm - 2 s — 1 . 

d The formal standard deviation for Fq in units of 
photonscm -2 s — 1 . 

e The redshifts z\ and z\§ at the values oi z at the mean 
fluxes for the first and tenth data bins from the left in 
Fig. 7 of 0.79 cm -2 s — 1 and 10.5 cm -2 s — 1 respectively. 
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TABLE 2 

MONOLUMINOSITY, q = 0.1 







Reduced x 2b 


F C 


<?F d 


l + ^i 

1 + Z10 


— u. 1U 


1U 


1 Q 


A 1 
4.1 


1 . 1 


1 OA 


-0.10 


12 


1.5 


4.0 


1.0 


2.00 


-0.10 


14 


1.2 


4.0 


1.0 


2.07 


-0.10 


16 


1.1 


3.9 


0.9 


2.15 


-0.10 


18 


1.1 


3.9 


0.9 


2.24 


-0.10 


20 


1.2 


3.8 


0.8 


2.34 


-0.10 


22 


1.4 


3.8 


0.7 


2.46 


-0.10 


24 


1.7 


3.8 


0.7 


2.60 


— U. 1 Z 


i n 
±u 


9 I 


a n 

■4:.U 


1 1 
± . ± 


1 87 


-0.12 


12 


1.9 


3.9 


1.1 


1.91 


-0.12 


14 


1.6 


3.8 


1.0 


1.95 


-0.12 


16 


1.3 


3.8 


1.0 


1.99 


-0.12 


18 


1.2 


3.8 


0.9 


2.05 


-0.12 


20 


1.1 


3.7 


0.9 


2.10 


-0.12 


22 


1.1 


3.7 


0.8 


2.16 


-0.12 


24 


1.1 


3.7 


0.8 


2.23 


n i a 

— U. 


i n 
±u 


9 7 


O.O 


1 1 
± . ± 


± .OU 


-0.14 


12 


2.4 


3.7 


1.1 


1.83 


-0.14 


14 


2.1 


3.7 


1.0 


1.85 


-0.14 


16 


1.8 


3.6 


1.0 


1.87 


-0.14 


18 


1.6 


3.6 


1.0 


1.90 


-0.14 


20 


1.4 


3.5 


0.9 


1.93 


-0.14 


22 


1.3 


3.5 


0.9 


1.96 


-0.14 


24 


1.2 


3.5 


0.9 


1.99 


— U. ID 


1U 


9 9 


9 A 


U.o 


i 7n 
1 . t u 


-0.16 


12 


2.1 


2.6 


0.8 


1.71 


-0.16 


14 


1.9 


2.6 


0.8 


1.72 


-0.16 


16 


1.8 


2.6 


0.8 


1.73 


-0.16 


18 


1.7 


2.6 


0.8 


1.74 


-0.16 


20 


1.5 


2.6 


0.8 


1.76 


-0.16 


22 


1.4 


2.6 


0.8 


1.77 


-0.16 


24 


1.3 


2.6 


0.8 


1.78 


-0.18 


10 


2.6 


2.4 


0.8 


1.65 


-0.18 


12 


2.5 


2.4 


0.8 


1.66 


-0.18 


14 


2.5 


2.4 


0.8 


1.66 


-0.18 


16 


2.4 


2.4 


0.8 


1.66 


-0.18 


18 


2.4 


2.4 


0.8 


1.66 


-0.18 


20 


2.3 


2.4 


0.8 


1.67 


-0.18 


22 


2.3 


2.4 


0.8 


1.67 


-0.18 


24 


2.2 


2.4 


0.8 


1.67 



a Spectral parameters that are defined in §3 and held 
constant in the x 2 minimization. 
b For 9 degrees of freedom. 

c The best fit value of the peak flux normaliza- 
tion in the 50keV to 300keV energy band in units of 
photonscm - 2 s — 1 . 

d The formal standard deviation for Fq in units of 
photonscm -2 s — 1 . 

e The redshifts z\ and z\§ at the values of z at the mean 
fluxes for the first and tenth data bins from the left in 
Fig. 7 of 0.79 cm -2 s — 1 and 10.5 cm -2 s — 1 respectively. 
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Fig. 7 — Examples of monoluminous peak-flux distributions fit to BATSE 
data. In the upper figures, the solid curve gives the BATSE data with 
statistical errors and the dotted curve gives the expected value from the best 
fit model. The peak flux is for the 256 ms timescale in the 50 keV to 300 keV 
energy band. The lower figures give the value of function of peak flux 

for the best fit model. For both models, \1/ = —0.14 and qo = 0.5. On the left, 
the reduced y 2 is 1.1 for 9 degrees of freedom with r T = 20, and, on the right, 
the reduced x 2 is 1-7 with t t = 10. 
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5.1. Monoluminosity Distribution 

In fits of the monoluminous flux distribution to the BATSE data, I fix the values of go, 
and tt- These variables are set by data that is independent of the peak flux distribution. 
In particular, the values of ^ and t t are set by fitting gamma-ray burst spectra. Generally 
one finds —0.16 < * < —0.12 and 10 < r T < 30 from fits to the spectra of bursts observed 
by BATSE (Brainerd et al. 1996b). This leaves two free parameters in the x 2 minimization, 
F and No, giving 9 degrees of freedom in the minimization. The results are given in Tables 
1 and 2 for \1/ ranging from —0.10 to —0.18 and Tt ranging from 10 to 26. In these tables, 
the first and second columns give the values at which ^ and are held, the third column 
gives the value of reduced x 2 from the minimization, the fourth column gives the value of 
Fq derived from the minimization, and the fifth column gives the formal error on Fq. The 
final column gives the ratio of 1 + z± to 1 + z w , where, counting from the lowest to the 
highest peak-flux rate, z\ is the redshift at the center of the first data bin and z w is the 
redshift at the center of the tenth data bin. These bins are centered on the peak-flux rates 
of 0.79cm~ 2 s _1 and 10.5 cm -2 s -1 respectively. The minimum values of reduced x 2 given in 
the tables are plotted against t t in Figure 6. 

The minimum value of x 2 is strongly dependent on the values chosen for \1> and r T , and 
only for specific combinations of these values are good fits found. The best value of reduced 
X 2 is w 1.0, which has an expectation of 0.5 for 9 degrees of freedom; an expectation of 
~ 0.01 is found for a reduced x 2 of 2.4. For a given value of q , the value of t t that 
gives a good fit increases rapidly as \1/ decrease, that is, as the value of E p decreases. 
The characteristic values of tt and ^ one derives from fits to the spectra (tt ~ 20 and 
\l/ m —0.14) produce a model fit to the flux distribution with a reduced x 2 °f 1-1 when 
go = 0.5, and a reduced x 2 °f 1-9 when g = 0.1. The values of F that one derives vary 
little as \1/ and r T vary; it is l-2cm _2 s _1 for g = 0.5, and it is 2.5-4 cm -2 s -1 for g = 0.1. 
As a consequence, the value of z at the BATSE threshold is 1.5-2.5. 

Figure 7 shows two specific examples of model fits, one of which fits the data well, 
the other poorly. When the model fits the data poorly, it is because the model breaks 
more rapidly away from a —5/2 power law than does the data. These figures also show the 
variation of z with F. 



5.2. Cutoff Power-Law Luminosity Distribution 

The flux distribution for the power-law luminosity distribution has three free 
parameters: F , N , and (5. Because go, \t r , and tt are fixed parameters, the model fits have 
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Fig. 8 — Reduced x 2 f° r fits of 
power-law luminosity peak-flux 
distributions to BATSE data. 
As in Fig. 6. 
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TABLE 3 
Power Law Luminosity, q = 0.5 



\]> a 


T T 


rvctlLlCctl \ 


-^0 


nr. d 


P 


!T P 


TVl f D f R 
Ullation JrxatlO° 




-0.10 


10 


0.9 


6.3 


9.9 


1.41 


0.70 


1.55 


1.38 


— 0.10 


12 


0.9 


6.3 


9.6 


1 .40 


0.77 


1 .62 


1 .43 


-0.10 


14 


1.0 


6.5 


9.1 


1.37 


0.83 


1.71 


1.48 


-0.10 


16 


1.2 


6.6 


8.2 


1.32 


0.86 


1.82 


1.55 


-0.10 


18 


1.3 


6.8 


6.9 


1.27 


0.84 


1.96 


1.63 


-0.10 


20 


1.6 


6.9 


5.5 


1.22 


0.77 


2.11 


1.72 


-0.10 


22 


1.9 


7.2 


4.4 


1.19 


0.64 


2.28 


1.80 


-0.10 


24 


2.2 


7.2 


3.4 


1.16 


0.54 


2.44 


1.88 


-0.12 


10 


0.8 


5.9 


10.0 


1.46 


0.63 


1.47 


1.33 


— 0.12 


12 


0.9 


5.9 


10.0 


1 .45 


0.68 


1.51 


1 .35 


-0.12 


14 


0.9 


6.0 


9.8 


1.43 


0.74 


1.56 


1.39 


-0.12 


16 


1.0 


6.2 


9.5 


1.41 


0.79 


1.62 


1.43 


-0.12 


18 


1.1 


6.3 


8.9 


1.37 


0.83 


1.69 


1.47 


-0.12 


20 


1.1 


6.5 


8.1 


1.34 


0.84 


1.77 


1.53 


-0.12 


22 


1.3 


6.6 


7.1 


1.30 


0.82 


1.87 


1.58 


-0.12 


24 


1.4 


6.8 


6.1 


1.26 


0.78 


1.97 


1.64 


-0.14 


10 


0.8 


5.5 


10.1 


1.50 


0.57 


1.40 


1.28 


— 0.14 


12 


0.8 


5.6 


10.1 


1 .49 


0.60 


1 .42 


1 .30 


-0.14 


14 


0.8 


5.6 


10.2 


1.49 


0.64 


1.45 


1.31 


-0.14 


16 


0.9 


5.7 


10.1 


1.47 


0.68 


1.48 


1.33 


-0.14 


18 


0.9 


5.8 


10.0 


1.46 


0.72 


1.51 


1.36 


-0.14 


20 


0.9 


5.9 


9.8 


1.44 


0.75 


1.55 


1.38 


-0.14 


22 


1.0 


6.0 


9.5 


1.41 


0.78 


1.60 


1.41 


-0.14 


24 


1.0 


6.2 


9.0 


1.39 


0.80 


1.64 


1.45 


-0.16 


10 


0.8 


5.2 


10.1 


1.53 


0.51 


1.35 


1.24 


-0.16 


12 


0.8 


5.2 


10.2 


1.53 


0.53 


1.36 


1.25 


-0.16 


14 


0.8 


5.3 


10.3 


1.53 


0.55 


1.37 


1.26 


-0.16 


16 


0.8 


5.3 


10.4 


1.52 


0.57 


1.38 


1.27 


-0.16 


18 


0.8 


5.4 


10.5 


1.52 


0.60 


1.39 


1.28 


-0.16 


20 


0.8 


5.4 


10.5 


1.51 


0.62 


1.41 


1.29 


-0.16 


22 


0.8 


5.5 


10.6 


1.50 


0.65 


1.43 


1.30 


-0.16 


24 


0.8 


5.5 


10.6 


1.49 


0.67 


1.45 


1.31 


-0.18 


10 


0.8 


5.2 


11.0 


1.58 


0.44 


1.29 


1.20 


-0.18 


12 


0.7 


5.2 


11.0 


1.58 


0.45 


1.29 


1.20 


-0.18 


14 


0.7 


5.3 


11.2 


1.58 


0.45 


1.30 


1.21 


-0.18 


16 


0.7 


5.3 


11.3 


1.58 


0.46 


1.30 


1.21 


-0.18 


18 


0.7 


5.3 


11.5 


1.58 


0.47 


1.30 


1.21 


-0.18 


20 


0.7 


5.3 


11.6 


1.58 


0.48 


1.31 


1.21 


-0.18 


22 


0.7 


5.3 


11.8 


1.58 


0.49 


1.31 


1.21 


-0.18 


24 


0.7 


5.3 


11.9 


1.57 


0.50 


1.31 


1.22 



a Spectral parameters that are defined in §3 and held constant in the \ 2 minimization. 
b For 8 degrees of freedom. 

c The best fit value of the peak flux normalization in the 50keV to 300keV energy band in units of 
photons cm -2 s — 1 . 

d The formal standard deviation for Fq in units of photons cm -2 s — 1 . 
e The best fit value of the luminosity function power law index. 
f The formal standard deviation for j3. 

s The ratio of the average time dilation for the mean fluxes of the first and tenth bins of Fig. 9 of 
0.79 cm - 2 s — 1 and 10.5 cm -2 s — 1 respectively. 

h The ratio of average spectral redshift for the mean fluxes used for the Dilation Ratio. 
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TABLE 4 
Power Law Luminosity, q = 0.1 



\]> a 


T T 


rvctlLlCctl \ 


-^0 


nr. d 


P 




Ullation JrxatlO° 




-0.10 


10 


0.8 


7.1 


15.2 


1.31 


1.08 


1.57 


1.44 


— 0.10 


12 


0.8 


7.2 


15.1 


1.29 


1.21 


1 .62 


1 .49 


-0.10 


14 


0.9 


7.3 


14.9 


1.24 


1.38 


1.70 


1.55 


-0.10 


16 


0.9 


8.1 


14.4 


1.08 


1.69 


1.85 


1.68 


-0.10 


18 


1.0 


8.6 


13.5 


0.96 


1.90 
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a Spectral parameters that are defined in §3 and held constant in the \ 2 minimization. 
b For 8 degrees of freedom. 

c The best fit value of the peak flux normalization in the 50keV to 300keV energy band in units of 
photons cm -2 s — 1 . 

d The formal standard deviation for Fq in units of photons cm -2 s — 1 . 
e The best fit value of the luminosity function power law index. 
f The formal standard deviation for j3. 

s The ratio of the average time dilation for the mean fluxes of the first and tenth bins of Fig. 9 of 
0.79 cm - 2 s — 1 and 10.5 cm -2 s — 1 respectively. 

h The ratio of average spectral redshift for the mean fluxes used for the Dilation Ratio. 
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Fig. 9 — Examples of power-law luminosity peak-flux distributions fit to 
BATSE data. Upper figures are as in Fig. 7. In the lower figures, the 
value of z m ax is given by *, the value of (z) is given by +, and the value 
of 1/ ((I + zy 1 ) - 1 is given by O. On the left, the reduced x 2 is 0.9 for 8 
degrees of freedom with t t = 20, and on the right, it is 0.8 with t t = 10. The 
value of z m in for the lowest bin in the lower left figure is 3.39, which is off the 
plot. 
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8 degrees of freedom. The introduction of a power-law luminosity distribution improves 
the ability of a cosmological gamma-ray burst model to fit the observed flux distribution, 
with most values of tt, and go producing acceptable values of reduced x 2 ■ These values 
are plotted in Figure 8, and the results of the model fits are listed in Tables 3 and 4. The 
first five columns in these tables are as in Tables 1 and 2. The sixth column is the derived 
value of f3, the seventh column is the formal error for f3, the eight column is the ratio of the 
average time dilation factor for the first data bin to the average time dilation factor for the 
tenth data bin, and the tenth column is the ratio of the average redshift factor for the first 
data bin to the average redshift factor for the tenth data bin. Two of the model fits listed 
in Table 3 are plotted in Figure 9. For the values of ^ and t t consistent with observed 
gamma-ray burst spectra, fits to the flux distribution produce values of reduced x 2 ~ 0.8. 

The derived values of F range between 5cm~ 2 s -1 and 10cm~ 2 s _1 , with the best 
values of x 2 having the lowest values of F . But the formal errors for F are of order 
10cm _2 s _1 . Such values for Fq imply isotropic burst luminosities of ~ 2 x 10 51 ergss~ 1 
when Hq = 75 kms" 1 Mpc -1 in the 50keV to 300 keV energy band. The value of (3 ranges 
between 0.7 and 1.6, with the highest values occurring for the lowest values of x 2 ■, but again 
the error is large, of order 0.5. These values of (3 are somewhat smaller than the value of 
1.9 one would get from a pure power-law fit to the data. The data therefore do not go 
to sufficiently low flux to exhibit the asymptotic limit. The values of (3 found in this study 
imply that the number of bursts per decade luminosity changes slowly with luminosity. 
For instance, with (3 = 1.6, only 75% of the bursts are within a factor of 10 of the lower 
threshold. This shows that gamma-ray bursts do not need a steeply falling or a narrow 
luminosity distribution to reproduce the observations. 

The large values of (3 in these fits imply that the redshift distribution function is 
strongly peaked at z < 1. The effect of this is clear in the lower halves of Figures 9a and 
9b. The value of z max for the bursts at threshold is generally ^ 3. On the other hand, 
the average time dilation and redshift factors are relatively small and slowly varying. In 
Figure 9, the implied redshifts from averages over (1 + z) and 1/(1 + 2;) are plotted. Even 
though z max ^ 3, the average redshifts are <^ 1. Comparing these figures to Figure 7 shows 
that the variation of the average z with F for the power-law luminosity model is less than 
the variation of z for the monoluminous model. For instance, in Figure 7a, the lowest flux 
bin, which is centered at 0.79 cm -2 s -1 , has z = 1.62 and the tenth bin, which is centered at 
10.5 cur 2 s^ 1 , has z = 0.37; in contrast, in Figure 9a, the average values of z from the time 
dilation are 1.18 and 0.41 for these flux bins, and the average values of z from the spectral 
redshift are 0.91 and 0.38. 

The values of z derived from fits to the peak flux distribution are lower but consistent 
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with the values of z derived from fits of the Compton attenuation spectral model to the 
spectra of bursts observed by BATSE. Spectral fits to relatively bright bursts (peak fluxes 
> 10cm _2 s~ 1 ) produce 0.5 < z < 4 (Brainerd et al. 1996b). Given the large uncertainty in 
Fq, those results do not contradict the results presented above. 

Three conclusions are evident. First, fits of the data to the monoluminous model are 
more sensitive to the spectral parameters than are fits to the power-law luminosity models. 
Second, the maximum redshifts that are consistent with the data are much larger when a 
luminosity function is present than when bursts are monoluminous. Third, the spectral 
redshift and time dilation from the cosmological expansion are weakened by the presence of 
a luminosity distribution. 



6. Redshift and Time Dilation Tests 

Researchers have sought evidence of a cosmological time dilation and spectral redshift 
in sets of gamma-ray bursts by comparing the characteristics of bursts at one flux to 
those at another. Such studies derive the ratio of the average time dilation (or redshift) 
at these two fluxes rather than the average value of z (Mitrofanov et al. 1994, 1996; 
Fenimore & Bloom 1995). When the bursts are monoluminous, this ratio goes to a 
constant as F — > oo. From equation (7), with the spectrum set to a power law oc u 5 , the 
limiting value of the redshift ratio (1 + z dirn ) / (1 + z brigh t) for the fluxes F dim < F bright 
is (F bright /F dim ) 1/{2 ~ S) when 1 < z < l/q and (F bright / F dim y 1/S when z > l/q . For a 
flux ratio of F bright / F dim = 10 and a characteristic power-law index of 5 = —3, the first 
relationship gives 2.15 and the second gives 1.58. The upper plot of Figure 10a gives 
the behavior of (1 + z dim ) / (1 + z bright ) as a function of F dim for F bright / F dim = 10. The 
asymptotic behavior for small flux is evident in this figure. More striking is the maximum 
in each curve at F xs 3 x 10~ 2 F . The maximum is generally <; 4, with the value of the 
peak decreasing with \1> and q . The values of this ratio for the model fits discussed above 
are given in Tables 1 and 2 for F dim = Fi = 0.79 cm -2 s™ 1 , the flux at the center of the first 
flux bin, and F bright = F w = 10.5 cm" 2 s" 1 , the flux at the center of the tenth flux bin. The 
value of z at the center of bin 1 is zi, and its value at the center of bin 10 is z±q. Note that 
this ratio is ^ 2 in the tables. Figure 10b gives the results for F brig h t / F dim = 5. Lowering 
the ratio of fluxes dramatically decreases the redshift ratio curves. 

Introducing a power-law luminosity function lowers the value of the average time 
dilation and redshift ratios from those found in the monoluminous model. In the lower 
plot of Figure 10a, one sees that the ratios of average time dilation are much smaller than 
in the monoluminosity case and that the peak value of the ratio is generally ^1.7. The 
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Fig. 10 — (1 + z dim ) / (1 + z bright ) versus peak flux. The ratios of time dilation 
factors for various cosmological models are plotted. The cosmological redshifts 
Zdim and Zbright are the redshifts associated with F^ m and Fbright, where, on the 
left, F bright = 10 F dim and, on the right, F bright = 5 F dim . The upper figures are 
for the monoluminous model, with the solid lines for q = 0.5 and the dotted 
lines for q = 0.1. The lower figures are for the power-law luminosity model 
with /3 — 1.6 and the remaining parameters as in the upper figures. For these 
figures, Zdim and Zbright are average values of redshift. For all curves tt = 18, 
and, from top to bottom in each figure, the curves have \1/ = —0.1, —0.14, and 
—0.18. The flux is normalized so that F = 1, and it is measured in the 50 keV 
to 300 keV energy band. 
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time dilation and redshift ratios for the first and tenth bins for the model fits to the data 
are given in Tables 3 and 4. The time dilation and redshift ratios for values of reduced x 2 
near 0.8 are <^ 1.6. These values can be as low as ~ 1.20 for some spectral parameters that 
produced good model fits to the flux distribution. This demonstrates that the cosmological 
effects are never large when a luminosity distribution determines the shape of the flux 
distribution, and that often the cosmological effects are slight, even though many bursts in 
the model have z ^ 1. Again, going to a flux ratio of Fb rig ht/ Fdim = 5 dramatically lowers 
these ratio curves. 

Many studies report a correlation of the gamma-ray burst timescale with the peak flux, 
finding that the dim bursts are on average longer than the bright bursts. Examples of the 
tests producing positive results are wavelet analysis, the total normalized counts test (Norris 
et al. 1994), and the fitted pulse width test (Davis et al. 1994). One test, the peak aligned 
test (Mitrofanov et al. 1993) produces conflicting results, with some researchers finding a 
correlation in the BATSE data (Norris et al. 1994), and others finding no correlation in 
either the BATSE data (Mitrofanov et al. 1996) or the data from the APEX experiment on 
Phobos-2 (Mitrofanov et al. 1993). This is somewhat disheartening, because clear evidence 
of a correlation would give the discussion that follows more resonance, but there is no 
consensus on the issue, so I'll address both sides of it. 

Those timescale studies that report a correlation generally find a dramatic effect. Norris 
et al. (1994) find that BATSE gamma-ray bursts with peak count rates above 18000s" 1 
are ~ 2.25 times as long as bursts with peak count rates between 1400 s" 1 and 2400 s" 1 . 
The actual time dilation ratio required to produce such an effect is much larger, because 
gamma-ray burst durations are photon-energy dependent. If the duration is proportional to 
^-0.23 (Mitrofanov et al. 1996), then a time dilation ratio of 2.8 is required. From Figure 10, 
one sees that a ratio of this size is impossible to achieve in the power-law luminosity model, 
and it is possible in the monoluminous model for only part of the available parameter space. 
As Fenimore & Bloom (1995) note, the monoluminous models produce a time dilation ratio 
of the required size at Zdi m ^ 6, which is inconsistent with the flux distribution. When the 
time dilation ratios for a factor of 10 difference in flux are calculated for the models in 
Tables 1-4, one finds that the monoluminous model produces ratios between 1.5 and 2 and 
the power-law luminosity model produces ratios between 1.2 and 1.5. 

To be more precise about this, let us use the specific model of a go — 0.1 universe and 
a ^ = —0.14 and = 20 spectral model. For a power law luminosity, a fit to the peak flux 
distribution gives F Q = 6.7 cm -2 s" 1 and (5 = 1.36, while for a monoluminous distribution, 
F = 3.5 cm" 2 s" 1 . 

In Davis et al. (1994), time dilation ratios of 1.9l^, 1.8±g;||, and 1.6±g;l? are found 
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from three different statistical tests when comparing bursts with count rates between 
18000 and 250000 counts per second to those between 2400 and 4500 counts per second. 
These tests are all based on fitting with width of pulse structure within a burst, and are 
therefore dependent on the photon energy. Taking this lower count rate to be a photon 
flux of 0.45 cm 2 s -1 , and taking the lower edge of the upper interval to be 5.8 cm 2 s -1 , and 
assuming that the width of the pulses go as z/~ 23 , one finds dilation ratios of 1.52 for the 
monoluminous model and 1.31 for the power-law luminosity model. Comparing the same 
upper count band to a lower range of 1400 to 2400 counts per second, the tests give ratios of 
2-0i[u7> 2.2l[J;44, and 1.8lo!5i, while the monoluminous model gives 1.52, and the power-law 
luminosity model gives 1.41. These various tests produces results that are all larger than 
the monoluminous model by about one standard deviation, and larger than the power-law 
luminosity model by about one and one-half standard deviations. For the quoted errors, the 
monoluminous model produces a value that is one and one-half standard deviations above 
unity, and the power-law luminosity model produces a value that is less than one standard 
deviation from unity. 

A test that compares the time between peaks, which is not energy dependent to the 
extent that the width of a single peak is, finds dilation ratios of 2.18 and 2.20 on the 256 ms 
timescale for two different selection criterion, with significances compared to unity of 0.013 
and 0.0016 respectively (Norris et al. 1996). Assuming gaussian statistics, these significances 
correspond to 2.5 and 3.2 standard deviations from unity, so the standard deviations for 
these measures are 0.47 and 0.38 respectively. In the absence of energy dependence of the 
measured duration, the monoluminous model gives 1.73 and the power-law luminosity model 
gives 1.42. The observed result is of order a standard deviation above the monoluminous 
model and two standard deviations above the power-law luminosity model. 

The one test that finds no evidence of a duration-peak flux correlation is the analysis 
of Mitrofanov et al. (1996). The peak-aligned method is used in this article to compare 
the average duration of bursts with different peak fluxes. It tests a different aspect of 
burst timescale than those just discussed, so the absence of a signal does not necessarily 
contradict these other tests. There is a contradiction with Norris et al. (1994), which 
Mitrofanov et al. show to be a consequence of the selection criterion use by Norris et al.. 
Mitrofanov et al. report on several measures of burst duration. In the first, bursts with peak 
fluxes < 1 cm" 2 s" 1 are compared to bursts with peak flux > 1 cm -2 s -1 . The total widths 
of the averaged light profiles are 6.57 ± 0.09 s and 6.64 ± 0.1 s respectively. This implies a 
duration ratio of 0.99 ± 0.13. In comparison, one finds a time-dilation ratio of 1.22 for a 
monoluminous model with a width energy dependence of z/~ 23 , and a time-dilation ratio of 
1.12 for a power-law luminosity model with the same energy dependence. The data is two 
standard deviations from the first, but it is only one standard deviation from the second. In 
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the second, bursts with peak fluxes < 0.45 cur 2 s -1 are compared to bursts with peak flux 
> 2.5 cm -2 s _1 . These bursts give total widths of 6.50 ± 0.13 s and 6.21 ± 0.15 s respectively. 
The duration ratio is then 1.05 ± 0.20. Taking 0.45 cm" 2 s" 1 for the low flux, one finds 1.49 
for the monoluminous models and 1.27 for the power-law luminosity model. This data is a 
little over two standard deviations from the first model value and one standard deviation 
from the second model value. From this, one sees that the data can be considered consistent 
with a cosmological time dilation if bursts have a broad luminosity distribution function. 

One study attempted to find a correlation between the peak flux and E peak , the 
observed photon energy at which the vF v curve has a maximum (Mallozzi et al. 1995). The 
advantages of using E peak rather than burst timescale are that E peak is strictly proportional 
to 1/ (1 + z), and its value varies by less than a factor of 10, whereas the duration is energy 
dependent, and it varies in value by several orders of magnitude. The study found that 
the average value of E pea k increases with peak flux, which one expects from cosmological 
effects. But the implied redshift factor is large. The subset of BATSE gamma-ray bursts 
with T90 > 2 s produces a ratio of average E peak of 2.45t[J;gg for bursts with peak fluxes 
between 5.90 cm -2 s" 1 and 105.0 cm -2 s _1 relative to bursts with peak fluxes between 
0.95cm _2 s _1 and 1.30 cm -2 s -1 . In contrast, model calculations for the parameters given 
above produce 1.61 for the monoluminous model and 1.30 for the power law luminosity 
model, which are 1.3 and 1.7 standard deviations from the data. For a lower peak flux 
range of 0.95cm _2 s _1 to 1.30cm _2 s _1 , the ratio is 1.99±o;3g, while the models give 1.47 
and 1.23 for the monoluminous and the power law luminosity cases, which are 1.4 and 2.0 
standard deviations from the data. 

Taken individually, the tests that show a correlation of burst duration or characteristic 
spectral energy with peak flux are consistent with the models, but taken together, they are 
inconsistent, because all of these tests are systematically above the monoluminous model 
by one standard deviation and above the power-law luminosity model by two standard 
deviations. The reported correlations are therefore inconsistent with the theoretical models 
presented in this article. In conclusion, the reported correlations must be a consequence of 
an intrinsic correlation of E peak and burst duration with peak flux. No cosmological model 
can simultaneously fit the observed flux distribution and the claimed time dilation and 
redshift ratios, and no cosmological model with a power-law luminosity distribution can fit 
the claimed time dilation and redshift ratios for any value of F . 
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7. The Flux Distribution for Redshift-Selected Subsets 



The observations require an intrinsic correlation of E veak and time dilation with 
luminosity. The question is whether a direct measure of this correlation exists. The recent 
studies of the flux distribution of gamma-ray bursts selected by hardness ratio presents such 
an opportunity. 

The discussion that follows depends on the two critical assumptions underlying 
equation (8). First, for z <C 1, n (r) is a constant. Second, the gamma-ray burst 
characteristics must be independent of r, particularly for z < 1. These assumptions are 
reasonable given that the observed flux distribution goes to a power-law of index —5/2 for 
large flux, which is a natural consequence of these assumptions. 

The discussion in §3 and §4 shows that the power-law index of the flux distribution 
deviates from —5/2 at F . This flux is dependent on the spectral parameters through 
T (L,p) in the L argument of $ (L,p), but this dependence is weak, as demonstrated by 
the values for F derived from model fits and given in Tables 3 and 4. If there is a strong 
dependence of F on spectral parameter, it can only enter through the direct dependence of 
the distribution function <3> (L,p) on the spectral parameter p. In this case, L Q is a function 
of p, which makes F a function of p through equation (6). 

The dependence of F on L is graphically demonstrated in Figure 11, where E peak is 
plotted as a function of flux for specific values of L and E p , where E p is the value of E pea k 
when z — 0. The solid curves give the variation of E peak with F when E p decreases with 
L, the dotted curves give E pea k for different values of L and the same E p , and the dashed 
curves give E pea k for different values of E p and the same L. In these curves, 2 C 1 when 
E pea k is constant with F, so the flux distribution for a given E pea k is a power law of index 
—5/2. Regarding these sets of curves as three models of the dependence of L on E p , one 
sees that when E p decreases with L , the value of F for bursts selected by E peak decreases 
with E p . The solid curves exhibit a second property not yet touched upon — the crossing of 
curves. This means that bursts with a given F and E pea k can have a variety of values for L 
and z. 

For the spectral model in §2, the selection of bursts by E peak is equivalent to turning 
the spectral parameter \1> into a function of E p and z. Equation (8) becomes 



The function * (E p ) is given by * (E p ) = a' c (E p /mc 2 ) Ep/armc 2 , where a c is the 



dN 






dV(Ep) 



(22) 
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Fig. 11 — Models of E peak 
versus F for constant L. From 
top to bottom, the solid lines 
are for * = -0.12 and F = 1, 
= -0.14 and F = 0.1, and 
* = -0.16 and F = 0.01. 
From top to bottom, the dotted 
lines are for Fq = 0.01 and Fq = 
0.1 with ^ = -0.12. From top 
to bottom, the dashed lines are 
for \I> = -0.14 and -0.16 with 
F = 1 and F = 0.1. The 
flux is measured in the 50keV 
to 300 keV energy band. 

Flux 

Klein-Nishina cross section and the prime represents the first derivative with respect to 
E p /mc 2 . For E p /mc 2 > 0.1972, ^ (E p ) « [log (2E p /mc 2 ) - l/2]mc 2 /E p , so the first 
derivative of * (E p ) is (E p ) /dE p = [\og(2E p /mc 2 ) - 3/2}mc 2 /E 2 . Therefore, the 
derivative term does not dramatically change the dependence of equation (22) on z from 
that in equation (8). The important difference from equation (8) is the dependence of $ on 
Epeak (1 + z) in equation (22). The observed distribution of E peak in the BATSE catalog 
falls rapidly to zero for 300 keV < E peak < 1 MeV (Mallozzi et al. 1995). Because the E peak 
distribution is narrow, $ {L,p) must fall rapidly as E p goes to infinity. This means that 
an additional upper limit on the value of z is set by the value of E peak (1 + z) at which 
the integrand no longer contributes significantly to the integral. When the integrand peak 
occurs at Zq <C 1, the maximum value of z is set by the dependence of <£> on L, but when 
the peak occurs at zq 1, the maximum value of z is set by the dependence on E p . In the 
second case, the upper limit on z is independent of F. This limit sets an upper limit on L 
through the ratio F/J r (z,p); as F decreases, the upper limit on L decreases. 

These results explain the results of recent studies (Pizzichini 1994; Belli 1995; 
Kouveliotou et al. 1996; Pendleton et al. 1997) of the BATSE flux distribution for 
gamma-ray bursts selected by hardness ratio. The hardness ratio is a measure of the 
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spectral shape, and it is dependent on the cosmological redshift because of the curvature 
in gamma-ray burst spectra below 1 MeV. It is therefore a proxy measure of E, pea ^. These 
studies show that the soft gamma-ray bursts have a flux distribution that is a power law of 
index —5/2. The hard gamma-ray bursts, however, follow a curve that falls below a —5/2 
power law. This is the behavior expected when E p rises with luminosity. 

The study of the burst distribution on the E^e^-log F plane may yield the luminosity 
distribution of gamma-ray bursts as a function of E p . This would enable one to separate 
the intrinsic correlation of E p with luminosity from the cosmological redshift. The way one 
approaches this problem is to select the group of burst with the highest value of E pea k- As 
one goes to lower flux for such bursts, the brighter bursts redshifted out of the E peak range 
being examined, and no bursts enter this band through a redshift from a higher value of 
E pea k- As a consequence, the peak-flux distribution deviates from the —5/2 power law in a 
manner determined by the luminosity distribution function for E p = E pea k, and a luminosity 
distribution function for this value of E p can be derived. One can then subtract from the 
density of bursts with smaller values of E pea k the density of bursts with large values of E p 
that are strongly redshifted. One then selects bursts with a lower value of E pea k to repeat 
the process. The uncertainty in the feasibility of this analysis is in whether the BATSE 
data set is of sufficient size and quality to perform such a test. 



8. Summary of Major Points 

A monoluminous burst distribution produces good fits to the BATSE peak-flux 
distribution for very specific spectral parameters. The parameters that produce the best 
values of \ 2 are themselves strongly dependent on the value of go- The spectra that are 
nearly power laws do not produce good fits: a curved spectrum is required to reproduce 
the observations. The spectral parameters derived from fits to gamma-ray burst spectra 
produce peak-flux distributions that fit the observed distribution well. From model fits to 
the peak-flux distribution, one finds that gamma-ray bursts at z — 1 have peak fluxes of 
~ 1-2 cm" 2 s" 1 for g = 0.5 and 2.5-4 cm -2 s -1 for g = 0.1. I confirm the conclusion 
that the go — 0.5 cosmology cannot simultaneously reproduce the observed flux distribution 
and the reported time dilation-burst flux correlations (Fenimore & Bloom 1995), and I find 
that the g <C 0.5 cosmology can never reproduce the observed time dilation-burst flux 
correlation, independent of the fit to the peak flux distribution. 

When gamma-ray bursts have a distribution of luminosities, the shape of the low end 
of the peak-flux distribution can be determined by the shape of the luminosity distribution 
(Meszaros & Meszaros 1995). If the luminosity distribution below some value is proportional 
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to L^ 13 , then the asymptotic flux distribution as F — > is oc F^ 13 when (3 ^ 0.6 for q = 0.5 
and (3 ^ 1.15 for q — 0.1 when there is only mild source density evolution. The precise 
values on the right side of these inequalities are dependent on the burst spectrum and 
the density evolution of burst sources. When the inequalities are violated, the peak-flux 
distribution is that found for monoluminous models. Equations (12) and (16) give these 
limits and their dependence on source density evolution. The luminosity distribution can 
formally diverge as the luminosity goes to 0; below a critical luminosity, which is set by 
the peak-flux threshold of the observing instrument, the luminosity distribution no longer 
contributes to the peak-flux distribution. A power-law luminosity distribution oc L^ 13 with 
an upper cutoff at L produces peak-flux distributions that fit the BATSE 3B data well 
for most spectral parameters. One finds that (3 ~ 1.5. This shows that the luminosity 
distribution need not be a standard candle or a steep power law to produce a good fit to 
the data. 

When the power-law luminosity distribution defines the shape of the flux distribution, 
so that (3 ~ 1.5, the gamma-ray bursts at a given peak flux have a distribution of redshifts 
that is strongly peaked at z < 1. As (3 becomes smaller, the peak in the z distribution 
goes to infinity, and the distribution rises monotonically. The distribution goes to zero 
at some maximum value for z that is set by the peak flux. When (3 ~ 1.5, the average 
redshift is a weak function of peak flux when the peak flux is small. At the peak-flux 
threshold of BATSE, the maximum redshift that a burst can have is ^ 3, while the 
average redshift at threshold is < 1. From the model fits to the peak flux distribution, 
the distribution in z has an upper limit of z — 1 when F ^ 5cm" 2 s _1 , although this 
limit on Fq is poorly constrained. This implies that the most luminous gamma-ray bursts 
have an isotropic luminosity of ^ 2 x 10 51 ergss _1 in the 50keV to 300 keV energy band 
for H = 75kms _1 Mpc -1 . For a burst duration of 10 s, the fluence is ^ 2 x 10 52 ergss _1 , 
making solar mass objects implausible sources of gamma-ray bursts. 

A distinction in the literature is often made between narrow and broad luminosity 
distributions. By narrow, many authors mean either that the bursts are nearly 
monoluminous, which is a statement about the burst physics, or that a high percentage 
of observed bursts have nearly the same flux, which is a statement about the detector 
threshold. A more useful distinction is between luminosity-dominated flux distributions 
and density-dominated flux distributions. In the former, the shape of the flux distribution 
is determined by the shape of the luminosity distribution, and the average redshift, which is 
weakly dependent on peak flux, is ~ 1. In the latter, the shape of the peak-flux distribution 
is determined by the spatial distribution of gamma-ray bursts, and there is nearly a 
one-to-one correspondence between redshift and peak flux. 
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The large errors in the model parameters do not strongly constrain the distribution 
in z at a given peak flux. Because of this, large values of z found from model fits of the 
Compton attenuation spectrum to BATSE gamma-ray bursts are consistent with the values 
of z possible in the flux distribution produced by a cosmological model with a luminosity 
distribution function. 

Because one averages over the distribution in z when calculating the average burst 
duration and the characteristic photon energy at a given peak flux, the average time dilation 
factor does not equal the average redshift factor; the former is larger than the latter. From 
the model fits to the observed peak-flux distribution, one expects bursts with two different 
values of peak flux, one 10 times the other, to have ratios of average time dilations and ratios 
of average spectral redshifts of < 1.6. Independent of the best fit to the flux distribution, 
the maximum possible time dilation and redshift ratios are < 1.7. The reported correlations 
of Ep ea k and burst time scale with peak flux are much larger than this by 1.5 to 2 standard 
deviations, so they are not consequences of cosmological expansion in the models examined 
above. Within this model, only the addition of an intrinsic correlation with luminosity can 
create the observed correlations. An alternative not explored in this article is strong source 
evolution, with a ~ 3 (Reichart & Meszaros 1997). From equations (12) and (16), the 
flux distribution in this case is density-dominated. On the other hand, since cosmological 
expansion produces a modest effect when a luminosity distribution determines the peak 
flux distribution, the absence of a correlation cannot be take as evidence that gamma-ray 
bursts are not cosmological. 

The presence of an intrinsic correlation of E peak with flux may be tested by examining 
the flux distribution for bursts with different values of E peak , the energy at which the vF v 
curve has a maximum. If the intrinsic value of E peak decreases with luminosity, then the 
bursts with the smallest observed values of E peak will have a peak-flux distribution that falls 
away from a —5/2 power law at a lower value of F than do the bursts with a high value of 
Ep e ak- This is observed when gamma-ray bursts are selected by hardness ratio, which is a 
proxy for E peak . 

One justification given for more sensitive gamma-ray burst detectors is the measurement 
of the gamma-ray burst flux distribution at fluxes below the thresholds of current gamma-ray 
experiments. If gamma-ray bursts are at cosmological redshifts of z ~ 1, then such an effort 
is of limited value, because the behavior of the distribution is strongly model dependent, 
and is easily modeled by a luminosity distribution function. Of more value are large 
numbers of gamma-ray bursts with well characterized spectra. Such burst catalogs allow 
one to discern the dependence of burst luminosity on spectral parameters, which then 
allows one to disentangle the intrinsic correlation of burst luminosity with E peak from the 
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correlation from the cosmological expansion. It is this that justifies accumulating as large a 
data set as possible, and underlies the value of continued observations by the BATSE. 
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